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Analog  sonobuoy  records  shot  over  0-10  m.y.  old  crust  along  the  East  Pacific 
Rise  have  been  digitally  scanned,  scaled,  and  analyzed  for  upper  crustal  velocity  structure. 
Previous  analyses  of  the  original  analog  records  [e.g.,  Houtz  and  Ewing,  1976]  suggested 
that  the  average  velocity  of  layer  2 A  at  the  EPR  is  3.4  km/s  and  increases  with  age  to  about 
4.5  km/s  in  40-m.y.-old  crust.  Later  work  [Diebold  and  Carlson,  1993]  suggests  that 
arrivals  interpreted  by  Houtz  and  Ewing  [1976]  as  representing  layer  2A  headwaves  are 
actually  diffractions  and/or  converted  shear-waves.  By  digitally  processing  and  analyzing 
many  of  the  original  analog  records  shot  near  the  EPR,  we  have  shown  that  (1)  true  layer 
2A  headwave  arrivals  are  not  present  in  these  records,  (2)  the  layer  2A  retrograde  branch 
reported  in  ESPs  and  other  surveys  cannot  be  resolved,  and  (3)  layer  2B/3  velocities  and 
average  velocity  gradients  directly  determined  from  the  data  range  from  4.4  to  7.6  km/s  and 
1.1  to  2.2  Is,  respectively.  Because  direct  estimates  of  'layer  2A'  velocities  are  not 
possible,  the  linear  gradient  approximation  for  the  hidden  layer  [Ewing  and  Purdy,  1982] 
was  applied;  estimated  velocities  at  the  top  of  the  basement  section  range  from  2.0  to  3.1 
km/s  and  velocity  gradients  are  in  the  range  1.9  to  3.4  Is.  These  results  show  considerable 
scatter,  and  no  relationship  between  velocities  in  the  uppermost  part  of  the  oceanic  crust 
and  age  of  the  crust  can  be  resolved. 

Introduction 

During  the  1960s,  ships  of  the  Lamont  Doherty  Geological  Observatory  collected 
over  500,000  km  of  variable-offset  airgun  sonobuoy  profiles  from  the  world’s  oceans. 
This  dataset  is  the  most  complete  suite  of  data  characterizing  the  in-situ  velocity  structure  of 
the  upper  oceanic  crust  and  sediment  cover.  The  data-collection  technique  involved 
launching  a  passive  expendable  sonobuoy  while  the  ship  was  underway.  This  sonobuoy 
would  remain  relatively  stationary  in  the  water.  While  moving  at  a  speed  of  about  6  knots 
(11  km/hr),  the  ship  would  tow  a  25  cu.  in.  airgun  firing  at  11  second  intervals.  The 
seismic  signal  would  be  received  by  the  launched  sonobuoy  and  transmitted  back  to  the 
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ship  via  an  FM  transmitter-receiver  system  [LePichon  et  al,  1968].  In  addition,  a  second 
sonobuoy  was  towed  200  m  behind  the  air  gun  and  recorded  a  normal  incidence  profile. 
The  recorded  analog  data  was  bandpass  filtered  (5-20  Hz)  and  plotted  as  sign-bit  traces  on 
thermally  sensitive  paper.  Figure  1  shows  an  example  of  the  final  data  output. 

The  upper  oceanic  crust  consists  of  "layers"  2A  and  2B/2C  .  "Layers"  2 A  and 
2B/2C  are  generally  accepted  to  be  composed  of  pillow  basalts  and  sheeted  dikes, 
respectively.  The  concept  of  a  layered  crust  originated  from  ophiolite  studies  and  analyses 
of  early  refraction  surveys  [e.g.,  Raitt,  1963;  Talwani  et  al.,  1971;  Houtz  and  Ewing, 
1976;  Houtz,  1976;  Salisbury  and  Christensen,  1978].  With  the  advances  of  seismic 
acquisition  techniques  (e.g.,  multi-channel  seismic  methods  and  ocean-bottom  sources 
and/or  receivers)  and  analyses  of  the  resultant  data,  this  layered  crust  concept  has  been 
replaced  by  a  crust  whose  physical  properties  vary  continuously  with  depth  [Vera  et  al., 
1990-,  Rohr  et  al.,  1990;  Carlson  and  Herrick,  \990‘,  Vera  and  Diebold,  1994;  Johnson  and 
Semyan,  1994;  Kappus  et  al.,  1995]. 

There  have  been  many  published  interpretations  of  portions  of  this  sonobuoy  data 
[e.g.,  LePichon  et  al.,  196S-,  Houtz  et  al.,  1968,  1970;  Ewing  and  Purdy,  1982;  Vera  and 
Mutter,  1988]  but  only  two  published  interpretations  of  a  large  portion  of  the  existing  data 
[Houtz  et  al.,  1970;  Houtz  and  Ewing,  1976].  Houtz  and  Ewing  [1976]  interpreted 
approximately  220  paper  sections  using  the  slope-intercept  method  and  concluded  that  they 
observed  definite  upper  crustal  velocity  variations  as  a  function  of  plate  age;  namely,  that 
seismic  layer  2A  velocities  increase  “from  about  3.3  kih/s  at  the  ridge  crests  to  that  of  layer 
2B  on  crust  about  40  m.y.  old.”  Since  their  publication,  Houtz  and  Ewing’s  [1976]  results 
have  been  debated  [e.g.  Jacobson,  1992;  Diebold  and  Carlson,  1993;  Carlson  and 
Jacobson,  1994].  The  purpose  of  this  research  is  to  reinterpret  a  portion  of  the  sonobuoy 
dataset  and  to  reevaluate  the  results.  The  sonobuoy  records  used  in  this  study  are  from 
seismic  profiles  shot  over  0-10  M.y.  old  crust  from  the  plates  surrounding  the  East  Pacific 
Rise  (EPR)  between  20°S  and  20°N  (see  Figure  2). 
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Data 

From  the  sonobuoy  record  shown  in  Figure  1,  it  can  be  seen  that  the  data  is  quite 
noisy.  This  noisiness  can  be  attributed  to  the  weak  signal  sourced  by  a  small  25  cu.  in. 
airgun  and  received  at  the  sonobuoy.  During  the  course  of  data  recording,  the  gain  setting 
was  changed  as  an  attempt  to  increase  the  recording  of  the  weak  received  signals.  This 
gain  change  has  had  the  effect  of  increasing  not  only  the  signal,  but  also  the  noise.  This 
gain  change  can  be  seen  in  Figure  1  at  approximately  3.7  km  offset. 

The  analog  sonobuoy  seismograms  were  recorded  as  sign-bit  traces  (i.e.,  without 
amplitude).  The  only  information  contained  in  the  traces  are  positive  and  negative  values 
recorded  as  black  and  white  as  a  function  of  offset  and  twt.  Therefore,  waveform  analyses 
could  not  be  made  on  this  dataset.  Another  obvious  feature  of  the  data  are  the 
reverberations  seen  from  the  received  signals.  The  reverberations  were  deliberately 
engineered  and  produced  at  the  source  to  enhance  the  appearance  and  interpretability  of  the 
recorded  signals.  This  has  been  useful  in  picking  the  arrivals  on  the  records. 

In  addition  to  the  recording  enhancements  and  type  of  signal  recorded  on  the 
seismograms,  various  sets  of  signals  can  be  seen  in  the  data.  These  sets  can  be  separated 
into  three  categories:  1)  seafloor  reflections,  2)  direct  wave  arrivals,  and  3)  refracted  wave 
arrivals.  Each  category  has  one  trait  in  common;  namely,  that  their  twt  increases  as  a 
function  of  offset.  The  seafloor  reflections  are  seen  as  a  hyperbola,  concave  downward 
(towards  increasing  twt)  and  emerging  from  the  zero  o^ffset  trace.  The  direct  wave  arrivals 
are  first  seen  at  zero  offset  and  zero  twt,  and  are  linear  with  an  inverse  slope  equal  to  the 
speed  of  sound  in  water.  Any  deviations  in  the  ship's  course  while  it  steamed  away  from 
the  sonobuoy  can  be  seen  as  departures  from  linearity  in  the  direct  wave  arrivals.  The  third 
category,  the  refracted  arrivals,  are  a  bit  more  subtle  on  the  sonobuoy  records.  They  are 
clearly  visible  only  as  branches  emanating  from  the  seafloor  reflection  and  are  concave 
upward  (i.e.,  toward  decreasing  twt).  These  refracted  arrivals  are  generally  not  tangent  to 
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the  seafloor  reflection  at  their  point  of  intersection.  This  implies  that  the  first  refracted 
wave  signal  that  was  received  and  is  observed  adjacent  to  the  seafloor  reflection  has  been 
turned  at  a  depth  significantly  greater  than  the  seafloor.  Thus,  waves  turned  at  shallower 
depths  would  be  seen  as  secondary  arrivals  behind,  or  later  than,  the  refracted  arrivals 
and/or  the  seafloor  reflections  [cf.  Stephen,  1982].  These  secondary  arrivals  from  shallow 
turning  waves  are  not  readily  visible  in  the  sonobuoy  records. 

Analysis 

Initially,  the  sonobuoy  records  were  digitally  scanned  and  scaled  at  Lamont- 
Doherty  Geological  Observatory.  These  digital  images  were  then  interpreted  within  a 
commercially  available  software  package.  Interpretation  of  the  images  on  the  computer 
involved  digitization  of  each  of  the  three  categories  of  the  seismic  data  mentioned  above 
(i.e.,  the  seafloor  reflections,  the  direct  wave  arrivals,  and  the  refracted  wave  arrivals). 
These  digitized  points  (offset-twt  pairs)  were  then  used  to  construct  oceanic  crustal  velocity 
models  based  on  the  assumption  of  stacked,  linear  velocity  gradient  layers.  These  velocity 
gradients  were  allowed  to  vary  as  a  function  of  depth.  In  addition,  the  model  layers  were 
assumed  to  parallel  the  seafloor  and  were  assumed  to  be  laterally  homogeneous.  These 
simplifying  assumptions  seem  reasonable,  considering  the  single  source-receiver 
configuration  of  the  data  collection  and  the  non-uniqueness  of  modeling  more  complicated 
models. 

Initially,  the  human  errors  in  picking  the  precise  locations  of  the  arrivals  of  the 
signals  on  the  sonobuoy  records  were  determined.  This  was  performed  by  picking  a  single 
observed  arrival  ten  times  at  the  same  offset  and  then  repeating  the  procedure  at  three  other 
offsets.  A  plot  of  the  digitized  picks  and  their  associated  standard  deviations  in  precision, 
or  repeatability,  of  picking  is  shown  in  Figure  3.  The  conclusion  here  is  that  the  average 
standard  deviation  in  the  human  picking  precision  is  slightly  less  than  3.0  ms.  With  this  in 
mind,  a  check  on  the  linearity  of  the  direct  wave  arrival  can  now  be  made.  This  linearity, 
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or  lack  thereof,  will  allow  for  appropriate  static  shifts  to  the  data  to  be  made  if  the  observed 
deviations  from  linearity  in  the  direct  wave  arrivals  are  greater  than  or  equal  to  the  error  in 
human  digital  picking.  At  some  offsets,  the  deviation  in  linearity  observed  for  the  direct 
wave  arrivals  is  greater  than  the  human  digital  picking  error.  This  means,  that  based  on  the 
analysis  of  the  direct  wave  arrivals,  appropriate  static  shifts  to  the  sonobuoy  data  can  be 
made.  Static  shifts  have  been  applied  to  individual  traces  to  correct  for  deviations  in  the 
course  of  the  ship  and  drift  of  the  sonobuoy. 

The  seafloor  reflection  is  the  next  dataset  category  to  be  analyzed  from  the 
sonobuoy  records.  The  seafloor  reflection  curve  was  used  to  model  the  normal  incidence 
water  depth  and  the  apparent  average  dip  of  the  seafloor  over  each  seismic  survey.  The 
average  absolute  dips  of  the  seafloor  range  from  0.0°  -5.1°  with  a  mean  of  1 .5°±1 .4°.  The 
average  dip  of  the  seafloor  calculated  from  each  sonobuoy  record  will  be  taken  into  account 
when  modeling  the  velocity  structure  of  the  crustal  layers. 

The  velocity  structure  of  the  crustal  layers  can  be  modeled  from  analyses  of  the 
refracted  wave  arrivals.  This  third  category  of  data,  the  refracted  arrivals,  contains 
information  about  the  velocity  structure  and  possible  layering  of  the  upper  oceanic  crust. 
The  basic  question  of  concern  here  is:  Assuming  parallel  layering  and  lateral  homogeneity, 
what  upper  crustal  velocity  structure  model  could  give  rise  to  the  observed  refracted 
arrivals?  As  mentioned  earlier,  there  is  generally  a  lack  of  tangency  between  the  reflected 
and  refracted  arrivals  at  their  point  of  intersection.  This  implies  that  waves  turned  just 
below  the  seafloor  are  not  being  observed  as  first  ‘arrivals  on  many  of  the  sonobuoy 
records.  If  this  is  the  case,  then  arrivals  from  these  shallow  turned  waves  may  be  second 
or  third  arrivals  and  thus  appear  to  be  masked  by  the  seafloor  reflection  and/or  the  primary 
refracted  arrivals.  Since  the  airgun  source  was  deliberately  engineered  to  provide  air 
bubble  oscillations  and  to  give  rise  to  the  apparent  ringing  in  the  data,  the  second  or  higher- 
order  arrivals  are  not  clearly  visible.  As  an  attempt  to  separate  the  first  arrivals  from  the 
higher-order  arrivals,  some  of  the  sonobuoy  records  were  transformed  to  the  tau-p  domain. 
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An  example  of  a  sonobuoy  record  transformed  to  the  tau-p  domain  is  shown  in  Figure  4 
[cf.  Vera  and  Mutter,  1988].  It  is  clear  from  this  figure  that  not  much  information,  if  any, 
has  been  gained  through  this  transformation  of  the  data.  The  noisiness  and  ringing  in  the 
data  appears  to  have  sufficiently  masked  these  higher-order  arrivals.  Thus,  arrivals  from 
waves  turned  at  shallow  depths  are  not  observed  on  the  sonobuoy  records.  To  model  the 
upper  crustal  velocity  structure  from  this  dataset,  assumptions  about  the  velocity  structure 
of  the  shallow  crust  where  arrivals  are  masked,  must  be  made.  Following  the  "hidden 
layer"  method  of  Ewing  and  Purdy  [1982],  it  is  assumed  that  the  uppermost  velocity 
structure  of  the  crust  consists  of  a  single  layer  with  a  linear  velocity  gradient  that  is  a 
function  of  depth. 

The  hidden  layer  approximation  is  a  simplistic  assumption  that  involves  modeling 
only  four  parameters,  the  thickness,  the  velocity  at  the  top,  the  velocity  at  the  bottom,  and 
the  linear  velocity  gradient  within  the  hidden  layer.  The  general  velocity  model  and  ray 
path  diagram  for  the  hidden  layer  model  is  shown  in  Figure  5.  A  straightforward 
calculation  of  these  parameters  is  discussed  in  Ewing  and  Purdy  [1982].  One  of  the 
parameters,  the  velocity  at  the  base  of  the  hidden  layer,  can  be  directly  measured  from  the 
sonobuoy  record.  It  is  equal  to  the  inverse  slope  of  the  refraction  branch  at  the  point  of 
intersection  of  the  reflected  and  refracted  arrivals.  The  other  three  parameters  can  be  found 
by  solving  the  following  set  of  equations: 


and 


T(X)-to  +  K;j  ln|^2vi  -  Kixi. 


"Ki  ■  Vl^Ki )  'y2) 
2hw  +  xosin0 

to  = — .  -  , 

VwV  1  -  P^Vw^ 

2hwpvw 

Xo=  - - = -  , 

yjl  -  p^Vw^  -  sin0 


(1) 

(2) 

(3) 

(4) 
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(modified  from  Nettleton  [1940])  where  p,  the  ray  parameter,  equals  the  slope  of  the 
refracted  branch  at  the  point  of  intersection  of  the  reflected  and  refracted  arrivals,  hw  is  the 
normal  incidence  water  depth,  is  the  velocity  in  water,  hi  is  the  normal  incidence 
thickness  of  the  hidden  layer,  X  and  T  are  the  offset  and  twt  of  the  point  of  intersection  of 
the  reflected  and  refracted  branches,  respectively,  while  xq  and  to  are  the  distance  traveled 
and  twt,  respectively,  for  the  portion  of  the  raypath  that  the  rays  traveled  in  the  water  layer, 
and  xi  is  the  distance  traveled  for  the  portion  of  the  raypath  that  the  wave  traveled  in  the 
hidden  layer  (first  layer).  0  is  the  dip  of  the  seafloor  in  radians  and  is  found  by  a  least- 
squares  fit  to  the  equation 


0  =  sin*^ 


Vw^t(Xr)^  +  4hw^ 
4hwXr 


(5) 


where  Xr  and  t(Xr)  are  the  offset  and  twt  vectors  for  points  along  the  seafloor  reflection 
curve. 

After  the  hidden  layer  velocity  structure  has  been  calculated,  the  deeper  or 
"observed"  layers  can  be  modeled.  These  layers  are  termed  "observed  layers"  because  the 
waves  that  were  turned  within  them  are  seen  as  first  arrivals  on  the  sonobuoy  records.  The 
number  of  observed  layers  giving  rise  to  the  refraction  data  is  difficult,  if  not  impossible,  to 
determine  from  the  data.  Therefore,  individually  stacked  layers,  each  having  a  linear 
velocity  gradient,  were  chosen  to  model  the  observed  layers.  The  number  of  layers  has 
been  determined  from  a  constrained  quadratic  spline  fit  to  the  data.  The  constraint  on  the 

i 

quadratic  spline  fit  is  that  the  second  derivative  is  less  than  or  equal  to  zero  at  all  points 
along  the  refracted  branch.  This  constraint  was  chosen  so  that  the  crustal  velocities  would 
never  be  allowed  to  decrease  with  depth  in  the  model;  i.e.,  low  velocity  zones  would  be 
prevented  from  being  present. 

The  number  of  observed  layers  was  chosen  from  the  number  of  knots  used  to  make 
the  constrained  quadratic  spline  fit  to  the  data.  There  were  generally  six  knots  needed  to 
adequately  fit  the  data,  while  attempting  not  to  fit  the  noise  in  the  data.  The  quadratic  spline 
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function  readily  lends  itself  to  calculating  the  linear  velocity  gradient  of  each  spline  piece. 
Six  spline  knots  gives  rise  to  five  spline  pieces,  or  five  observed  layers.  The  velocity 
gradient  (Kn),  velocity  at  the  top  (Vn.i),  velocity  at  the  bottom  (vn),  and  thickness  (hn)  of 
each  observed  layer  can  be  calculated  from  the  quadratic  function  for  each  spline  piece 
using  the  following  formulas: 


\i 

-  P^Vn-i^ 

/ 

X  =  xo  + 

where  xq  can  be  calculated  from  equation  (4)  and 


Xi  = 


2^1  -  p^Vj.i^ 


2vr^ 


2v.2 


p^v 


pKi 


(6) 


(7) 


where  Vj  =  Vj.i  +  Kihi. 

Also,  T(X)  =  to  4^tj  +  In 

i~\ 

where  to  can  be  calculated  from  equation  (3)  and 


1  +  Vl  -  p^Vn-l^ 
PVn 


(8) 


.(9) 


Equations  (6)  and  (8)  were  simultaneously  solved  for  Kn.  This  was  performed  in  the  tau-p 
domain  (i.e.,  tau  =  t  -  px)  using  the  iterative,  damped  least  squares  technique  of  Marquardt 
[1963].  The  smooth  spline  function  has  the  added  benefit  of  modeling  a  continuous 
velocity  function  throughout  the  layers.  Since  the  spline  is  a  smoothly  varying  function 
without  discontinuities,  then  the  velocity  field  (equal  to  the  inverse  slope  of  the  spline 
function)  is  also  without  discontinuities.  A  plot  of  the  digitized  refracted  branch  and  the 
constrained  quadratic  spline  fit  to  the  data  is  shown  in  Figure  6.  Also  shown,  is  the 
resultant  velocity  structure  with  associated  uncertainties  calculated  from  the  spline  fit  to  the 
data.  The  complete  list  of  sonobuoy  modeling  results  are  listed  in  the  Appendix. 
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Discussion 

Having  interpreted  over  200  sonobuoy  records  from  the  Pacific  and  Atlantic  using 
the  slope/intercept  method,  Houtz  and  Ewing  [1976]  have  proposed  that  (1)  the  seismic 
structure  of  layer  2  consists  of  three  sublayers  2A,  2B,  and  2C  having  p-wave  velocities  of 
3. 6+0.4,  5;2±0.4,  and  6.1±0.2  km/s,  respectively;  (2)  layer  2A  velocities  increase  with 
age  over  tens  of  millions  of  years  until  they  cannot  be  distinguished  from  velocities 
characteristic  of  layer  2B;  and  (3)  layer  2 A  apparently  thins  with  age,  disappearing  when 
velocities  in  layer  2A  approximate  those  in  layer  2B  after  about  40  m.y.  [Carlson  and 
Jacobson,  1994].  Their  three-layer  model  that  was  based  on  a  slope/intercept  interpretation 
has  some  similarities  with  the  model  and  analyses  that  has  been  demonstrated  here.  For 
example,  if  the  upper  crust  consisted  of  three  distinct  layers,  then  the  modeling  using  the 
constrained  quadratic  spline  outlined  here  would  have  shown  these  three  layers.  In 
general,  many  more  layers  are  need  to  model  the  upper  crust,  but  the  three  layer  model 
appears  to  fit  the  data  to  first  order  and  is  a  good  starting  model. 

The  suggestion  by  Houtz  and  Ewing  [1976]  that  layer  2A  (the  layer  at  the  top  of 
basement)  becomes  indistinguishable  from  layer  2B  at  approximately  40  m.y.  was  based 
for  the  most  part,  on  apparent  head  wave  arrivals  from  the  top  of  basement.  Diebold  and 
Carlson  [1993]  have  shown  that  many  of  these  interpreted  head  wave  arrivals  may  be 
diffraction  hyperbolas  or  converted  shear-wave  arrivals,  not  critically  refracted  head  waves 
from  the  top  of  layer  2A.  So,  the  age-dependent  variation  reported  by  Houtz  and  Ewing 
[1976]  may  not  be  entirely  accurate.  The  analyses  of  the  sonobuoy  records  presented  in 
this  report  can  be  checked  for  age-dependency  of  the  velocity  structure  of  the  uppermost 
oceanic  crust.  This  can  be  done  by  analyzing  the  age  variations  in  the  parameters 
determined  for  the  hidden  layer  in  this  report.  Graphs  of  the  four  hidden  layer  parameter 
(thickness,  velocity  at  the  top,  velocity  at  the  bottom,  and  linear  velocity  gradient)  versus 
age  are  shown  in  Figure  7.  Significant  age  variations  in  the  hidden  layer  parameters  do  not 
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appear  to  exist  from  the  analysis  of  the  total  sonobuoy  dataset  presented  here.  This  could 
be  explained  by  considering  that  oceanic  crust  formed  at  different  portions  along  the  EPR 
may  be  mineralogically  and  mechanically  different,  thus  overprinting  any  systematic  age 
dependent  velocity  variations.  Perhaps  a  better  check  on  age  variations  would  be  to  analyze 
the  velocity  modeling  results  from  a  single  transect  running  perpendicular  to  the  ridge 
crest.  This  type  of  transect  should  sample  crust  that  was  formed  at  relatively  the  same  place 
along  the  EPR.  Referring  back  to  the  sonobuoy  basemap  (Figure  2),  there  is  a  transect  that 
runs  perpendicular  to  the  ridge  crest.  This  transect  is  located  at  approximately  12°45'N 
latitude  and  sampled  crust  ranging  in  age  from  1  -  9  m.y.  Figure  8  shows  graphs  of  the 
hidden  layer  parameters  versus  age  from  this  transect.  Again,  there  does  not  appear  to  be 
any  statistical  age  variation  in  the  uppermost  crustal  velocity  stmcture. 

The  12“45’N  transect  velocity  modeling  results  can  also  be  used  to  check  on  age 
dependent  variations  in  the  deeper  (observed)  layers.  If  the  velocity  structure  changes  with 
age,  then  an  apparent  thinning  or  thickening  of  velocity  isopachs  with  age  should  occur. 
For  example,  if  layer  2A  thins  with  age  (as  Houtz  and  Ewing  [1976]  stated),  then  the 
isopach  from  the  seafloor  to  the  3.6  km/s  (or  similarly,  4.0  km/s)  iso-velocity  surface 
should  thin  with  crustal  age.  Similar  analyses  can  be  done  for  other  isopachs.  A  plot  of 
the  depths  to  the  4.0,  5.0,  and  6.0  km/s  iso- velocity  surfaces  from  analyses  of  the  northern 
sonobuoy  transect  is  shown  in  Figure  9.  The  iso-velocity  lines  do  not  appear  to  show  any 
statistical  age  variation;  i.e.,  there  is  not  any  statistically  recognizable  thinning  or  thickening 
of  the  velocity  isopachs.  Thus,  from  the  analyses  of  this  sonobuoy  dataset,  there  does  not 
appear  to  be  any  statistical  age  variation  in  the  velocity  structure  for  this  portion  of  the 
upper  oceanic  crust.  This  result  is  consistent  with  those  stated  by  Purdy  [1982]  for  airgun- 
ocean  bottom  hydrophone  surveys  at  12°N  near  the  EPR  and  Vera  and  Mutter  [1988]  for 
sonobuoy  and  expanding  spread  profiles  over  young  crust  between  11°  and  13°N  near  the 


EPR. 
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In  the  past,  researchers  have  discussed  age  dependent  mechanisms  that  could  give 
rise  to  expected  age  dependent  velocity  variations  (e.g.,  Houtz  and  Ewing,  1976;  Schreiber 
and  Fox,  1977;  Carlson  and  Herrick,  1990;  Rohr,  1994;  Grevemeyer  and  Weigel,  1996). 
To  account  for  the  modeling  results  shown  here,  a  balance  of  age  dependent  mechanisms 
that  could  give  rise  to  the  lack  of  p-wave  velocity  variations  must  be  present.  For  example, 
infilling  of-voids  and  cracks  in  the  layer  2A  pillow  basalts  over  time  will  lead  to  a 
significant  density  increase,  thus  lowering  the  seismic  velocity;  while  at  the  same  time,  this 
added  sediment  and  alteration  products  will  act  to  stiffen  the  overall  basalt,  thus  raising  the 
seismic  velocity.  The  overall  effect  could  be  to  insignificantly  change  p-wave  velocity. 
Another  example  stems  from  the  effect  that  cracks  play  in  determining  seismic  velocities.  It 
has  been  shown  that  low  aspect  ratio  cracks  (flat  cracks)  lower  p-wave  velocities  much 
more  significantly  than  high  aspect  ratio  cracks  (spherical  cracks)  [e.g.,  Walsh,  1965; 
Kuster  and  Toksoz,  1974]. One  would  expect  to  find  a  higher  concentration  of  high  aspect 
ratio  cracks  in  young  basaltic  crust  (e.g.,  in  pillows  and  lava  tubes)  than  in  older,  more 
consolidated  basaltic  crust.  This  would  have  the  effect  of  both  reducing  the  total  porosity 
with  age  (raising  velocity)  and  also  changing  the  distribution  of  crack  aspect  ratios  more 
towards  the  lower  end,  flat  cracks  (lowering  velocities)  The  net  effect  of  reducing  porosity 
and  changing  the  overall  crack  aspect  ratio  distribution  with  age  may  also  be  to 
insignificantly  alter  the  p-Owave  velocity  of  the  upper  oceanic  crust.  These  are  just  two 
possible  explanations  for  the  observed  lack  of  age  dependency  in  the  observed  seismic 
velocities.  Other  mechanisms  undoubtedly  also  affect  !the  balance  of  p-wave  velocity  with 
age. 

Conclusion 

Sonobuoy  records  taken  from  20°S  -  20°N  along  the  East  Pacific  Rise  over  crust 
ranging  in  age  from  0  to  10  m.y.  have  been  digitally  scanned,  scaled,  and  modeled  for 
upper  crustal  velocity  structures.  Modeling  of  the  velocity  structure  consisted  of  analyses 
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of  shipcourse  navigation  problems,  dip  moveout  corrections  from  the  seafloor  reflections, 
and  refracted  or  turned  ray  analyses.  We  have  shown  that  (1)  true  layer  2A  headwave 
arrivals  are  not  present  in  these  records,  (2)  the  layer  2A  retrograde  branch  reported  in 
ESPs  and  other  surveys  cannot  be  resolved,  and  (3)  layer  2B/3  velocities  and  average 
velocity  gradients  directly  determined  from  the  data  range  from  4.4  to  7.6  km/s  and  1 . 1  to 
2.2  /s,  respectively.  The  modeling  results  show  a  lack  of  statistical  age  variations  in  the 
upper  crustal  p-wave  velocity  structure  for  this  region  along  the  East  Pacific  Rise.  Several 
reasons  could  account  for  the  lack  of  age  variation,  including  density  and  crack  aspect  ratio 
distribution  changes  with  age  for  the  upper  oceanic  crust. 
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Appendix 


vi  is  the  velocity  at  the  base  of  the  i-th  layer,  | 

Ki  is  the  velocity  gradient  within  the  i-th  layer,  and 
hi  is  the  thickness  of  the  i-th  layer. 
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Figure  Captions 

Fig.  1.  Variable-offset  airgun  sonobuoy  profile  105V28-B  shot  over  6.26  m.y.  old  crust  at 
95°15'W,  5'’19'N.  The  data  has  been  bandpass  filtered  (5  to  20  Hz)  and  is  composed  of 
sign-bit  (Os  and  Is)  traces.  Note  the  change  in  gain  setting  at  approximately  3.7  km  offset, 
and  the  loss  of  clarity  below  the  seafloor  reflection  and  ringiness  of  the  data  due  to  air- 
bubble  source  oscillations. 

Fig.  2.  Locations  of  the  variable-offset  airgun  sonobuoy  profiles  used  in  this  research. 
The  entire  dataset  lies  near  the  East  Pacific  Rise  over  crust  ranging  in  age  from  0  to 
approximately  10  m.y. 

Fig.  3.  Plot  of  digitized  picks  at  four  different  offsets  along  the  observed  refracted  branch 
on  the  sonobuoy  records.  The  refracted  branch  was  picked  ten  times  at  each  of  the  four 
offsets  to  determine  the  repeatability,  or  precision  in  signal  picking.  The  uncertainty  in 
picking  was  determined  from  the  average  standard  deviation  of  the  picked  values.  The 
human  picking  precision  is  slightly  less  than  3.0  ms. 

Fig.  4.  Tau-p  transform  of  the  sonobuoy  record  105V28-B  shown  in  Figure  1.  This 
transformation  was  done  as  an  attempt  to  better  display  the  refracted  arrivals  from  turned 
rays.  This  transformation  of  the  data  was  unsuccessful  in  showing  more  information  in  the 
records  than  can  be  seen  in  the  original  traveltime  vs.  offset  displays.  The  noise  in  the  data 
appears  to  be  even  more  evident  in  the  tau-p  transformation  than  in  the  original  records. 
Thus,  all  analyses  in  this  report  were  done  on  the  original  traveltime  vs.  offset  records. 
Compare  this  figure  with  Figure  3a  ofVera  and  Mutter  [1988]. 

Fig.  5.  (a)  Ray  path  diagram  and  velocity  model  for  the  hidden  layer  approximation  of 
Ewing  and  Purdy  [1982],  It  is  assumed  that  there  is  a  single  hidden  layer  with  a  linear 
velocity  gradient  with  depth  present  in  the  uppermost  oceanic  crust  whose  arrivals  are 
masked  by  the  seafloor  reflection  and  are  thus  uninterpretable,  (b)  Forward  model  on  the 
right  of  the  velocity  structure  shpwn  on  the  left.  Note  the  triplication  of  the  refracted 
branch  being  overshadowed  by  the  dominant  seafloor  reflection.  From  the  example  in 
Figure  1,  it  is  apparent  that  the  arrivals  from  rays  turned  at  shallow  depths  are  not  visible  in 
the  data. 
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Fig.  6.  Constrained  quadratic  spline  fit  to  the  digitized  refracted  branch  observed  in  the 
original  sonobuoy  record  (SB  105V28-B)  is  shown  in  (a).  This  quadratic  spline  readily 
allows  for  the  calculation  of  a  six  layer  crustal  velocity  model  -  one  hidden  layer  and  five 
layers  modeled  from  the  spline  pieces.  Each  layer  is  modeled  with  a  linear  velocity 
gradient,  (b)  Resultant  velocity  model  with  associated  uncertainties  represented  by  a 
velocity  envelope.  Any  velocity  function  that  would  fall  within  the  envelope  adequately  fits 
the  data. 

Fig.  7.  Calculated  parameters  for  the  hidden  layer  as  a  function  of  crustal  age  for  the  entire 
dataset  shown  on  the  basemap  in  Figure  2.  No  statistical  correlation  was  found  between 
any  of  these  parameters  and  age. 

Fig.  8.  Calculated  parameters  for  the  hidden  layer  as  a  function  of  crustal  age  for  the 
northern  transect  at  approximately  12°45'N  latitude  shown  on  the  basemap  in  Figure  2.  No 
statistical  correlation  was  found  between  any  of  these  parameters  and  age. 

Fig.  9.  Depth  to  the  4.0,  5.0,  and  6.0  km/s  constant  velocity  surfaces  (iso-velocity 
surfaces)  modeled  from  the  northern  transect  dataset.  From  this  figure,  no  apparent 
thinning  or  thickening  of  velocity  isopachs  with  crustal  age  exists.  Layer  2A  does  not 
appear  to  thin,  nor  does  layer  2B/2C  appear  to  thicken  away  from  the  ridge  crest.  Upper 
crustal  p-wave  velocity  structure  does  not  appear  to  be  a  function  of  crustal  age  for  the 
airgun  sonobuoy  records  analyzed  in  this  report. 
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